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We consider approaches to cosmological parameter estimation in the inflationary cosmology, fo- 
cussing on the required accuracy of the initial power spectra. Parametrizing the spectra, for example 
by power-laws, is well suited to testing the inflationary paradigm but will only correctly estimate 
cosmological parameters if the parametrization is sufficiently accurate, and we investigate conditions 
under which this is achieved both for present data and for upcoming satellite data. If inflation is 
favoured, reliable estimation of its physical parameters requires an alternative approach adopting 
its detailed predictions. For slow-roll inflation, we investigate the accuracy of the predicted spectra 
at first and second order in the slow-roll expansion (presenting the complete second-order correc- 
tions for the tensors for the first time). We find that within the presently-allowed parameter space, 
there are regions where it will be necessary to include second-order corrections to reach the accu- 
racy requirements of MAP and Planck satellite data. We end by proposing a data analysis pipeline 
appropriate for testing inflation and for cosmological parameter estimation from high-precision data. 

PACS numbers: 98.80.Cq, 98.70.Vc astro-ph/0202094 



I. INTRODUCTION 

Recent cosmic microwave background (CMB) 
anisotropy results Q] , showing a multiple peak structure 
in the anisotropy power spectrum, lend powerful support 
to the inflationary cosmology as the origin of structure in 
the Universe. It is now widely expected that cumulative 
improvements in the CMB data will lead to progressively 
more accurate estimation of cosmological parameters, 
with projects funded so far culminating in the Planck 
satellite mission expected to report results around 2010. 

Given a set of data on structures in the Universe, such 
as the CMB power spectrum, it is necessary to simulta- 
neously fit both for the parameters describing the global 
cosmology (such as the matter budget and expansion 
rate) and those describing the so-called 'initial pertur- 
bations'; they cannot be considered separately. If the 
model for the initial perturbations is insufficiently ac- 
curate, or even worse completely wrong, the full power 
of the experiment to constrain cosmological parameters 
cannot be exploited. 

The inflationary cosmology is an attractive paradigm 
for the generation of the initial perturbations, but even 
there the situation can be very complicated in general. 
In particular, if there are multiple scalar fields the per- 
turbations can be a mixture of isocurvature and adia- 
batic, and may be non-gaussian. Such initial conditions 
may prove difficult or even impossible to parametrize, 
and if such an inflation model is correct it will be a ma- 
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jor obstacle to successful parameter estimation. However 
it remains a powerful working hypothesis that the sim- 
plest class of models, where inflation is driven by a single 
scalar field, is viable; this creates a framework within 
which the necessary calculations are reasonably simple, 
with the initial perturbations computed either approxi- 
mately analytically or exactly numerically. As yet, there 
is no indication from observations that we might need to 
go beyond this paradigm. 

The main goal of this article is to investigate differ- 
ent strategies that an observer can use to estimate the 
cosmological parameters, and to examine the extent to 
which it is necessary to adopt detailed inflationary pre- 
dictions. The spectrum of the fluctuations is assumed to 
be produced by an underlying inflationary model and is 
calculated exactly by means of numerical computations. 
Given this situation, we study how the data analysis can 
be performed in two different scenarios. The first scenario 
applies if one wants only to estimate cosmological param- 
eters, such as the baryon density and reionization optical 
depth, and does not care about the underlying inflation 
model beyond being confident that the description of the 
initial perturbations used is adequate. In this case, ob- 
servers typically use a power-law fit, see e.g. Ref. pj, and 
the first question is to test how accurate a power-law fit is 
to typical inflationary cosmologies. In particular, we wish 
to know if this kind of fit is accurate enough for present 
data, and whether it will also be accurate enough to an- 
alyze high-precision data like that to be provided by the 
Planck satellite. The second scenario, which makes more 
stringent requirements on theoretical accuracy is if one 
intends to estimate properties of the inflationary model. 
In this case, the slow-roll method can be used to calculate 
an approximate spectrum and we will study its accuracy. 
We also consider to which order in the slow-roll parame- 
ters the spectrum should be calculated in order to reach 
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the Planck precision. We propose an analysis pipeline 
for testing the consistency of single-field slow-roll infla- 
tion and estimating physical parameters of inflation, e.g. 
the energy scale of inflation. 



II. INFLATIONARY BASICS 

A single-field inflation model generates Gaussian spec- 
tra of purely adiabatic density perturbations (scalar per- 
turbations) and gravitational waves (tensor perturba- 
tions). We denote the dimensionless power spectra by 
Vizik), where 1Z is the intrinsic curvature perturbation 
on comoving hypersurfaces (identical with Bardeen's £ Q 
up to a sign), and Vh(k), h being the amplitude of gravi- 
tational waves. Scalar and tensor perturbations obey the 
equation of a parametric oscillator Q 
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where a prime denotes differentiation with respect to 
conformal time rj and k is the comoving wavenum- 
ber. This equation only requires the assumption of 
linear perturbation theory. The quantities /xs,t(??) are 
defined by fisiv) = 2zslZ and ^t{v) = where 
zg = ay2 — aa"/a' 2 and zt = a. The initial conditions 
for the mode functions fig,T are fixed by the assumption 
that the quantum fields are in the vacuum state when 
the mode k is subhorizon, 
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where rji is an arbitrary initial time at the beginning of 
inflation. The power spectra are calculated according to 
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Both power spectra can be derived from the inflaton po- 
tential V(<f)) and the initial conditions for the inflaton 
field 0, and hence are not independent. They can be 
obtained numerically by solving the appropriate mode 
equations wavenumber by wavenumber (see e.g. Ref. Q). 
In the following, they are denoted by 'Pnum(fc)- 
The tensor-to-scalar ratio 
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is of interest for testing the consistency of a given model 
of inflation. It has often been defined in terms of the 
microwave background quadrupole moments. This defi- 
nition has the disadvantage that it depends on the cosmo- 
logical parameters, especially the density of the cosmo- 
logical constant f^A 0,0 . In Ref. Q the ratio between V$, 
(where $ is the gauge-invariant Bardeen metric potential 
|8j) and Vk was used, which removes the dependence on 



Da- However $ does still depend on the dynamics of the 
Universe at the photon decoupling epoch and thus is not 
completely model independent (it depends mainly on the 
physical matter density u) m = fl m h 2 ). The advantage of 
1Z is that it is conserved on super-horizon scales once the 
decaying mode is negligible and provided only adiabatic 
perturbations are considered 0, EJ El ■ 

The spectral indices and their "running" are defined 
by the following expressions 
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For purposes of illustration, in this paper we use three 
qualitatively different inflationary models to mimic ide- 
alized measurements of the power spectra. The first is a 
chaotic inflation model with a quartic potential 



V(<f>) = \4>\ 

the second a false vacuum inflation potential 
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with fi = 1, which is inspired by the scenario of hy- 
brid inflation and the third a potential introduced 
in Ref. El 
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For each potential we need to know the scalar field 
value 4>* when observable perturbations were generated 
(i.e. when a given scale fc* was equal to the Hubble radius 
during inflation), corresponding roughly to 55 e-foldings 
from the end of inflation. The last two potentials provide 
no natural end to inflation, and we make an arbitrary 
choice for 4>.j, to be equal 0.3\/2mpi and — 0.3mpi respec- 
tively. The chaotic inflation model ends by violation of 
slow-roll and so we take (/>* ~ 4.2mpi. 

Fig. n shows the scalar and tensor power spectra 
'Pnum(fc) for these three models, obtained numerically by 
the method of Ref. The corresponding microwave 
anisotropics, obtained using a modified version of CAMB 
[l4|. are also shown. 1 We present the characteristic quan- 
tities of these spectra, evaluated at k* , in Tabic HJ Wc 
have chosen these three models because the chaotic model 
is an example in which tensor perturbations are relevant 
and it shows moderate negative tilt, the false vacuum 
model has a moderate positive tilt and the arctan model 
has both large tilt and running. 



1 A module to directly input the predictions of slow-roll in- 
flation to the CAMB program is available to download at 
www. astronomy . susx . ac . uk/ ~sleach/ inflation/ 
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FIG. 1: The top panel shows the power spectra of scalar (up- 
per lines) and tensor (lower lines) perturbations for our three 
models. The scalar spectra are normalized to Vn — 2 x 1CP 9 
at the scale fc* = 0.01/iMpc -1 , which approximately matches 
the COBE normalization. The bottom panel shows the 
corresponding Ce curves for a flat cosmological model with 
io b = 0.0200, Lu m = 0.1268 and u A = 0.2958 (implying 
h — 0.65), and reionization optical depth r = 0.05, with 
the upper lines again the scalar contribution and the tensors 
considerably subdominant. Only the sum of the two can be 
detected, though they contribute differently to polarization 
anisotropies. 



III. PARAMETER ESTIMATION IGNORING 
INFLATIONARY PREDICTIONS 

A. Parametrizing the spectra 

To estimate the cosmological parameters we need an 
adequate parametrization of scalar perturbations, and 
more sophisticated analyses informed by inflation also in- 
clude tensors. 2 If one is only interested in a measurement 



Exact values R ns — 1 ht as ot 



Chaotic 


0.285 


-0.055 


-0.037 


-0.0009 


-0.0006 


False vacuum 


0.051 


0.054 


-0.006 


0.0018 


0.0005 


Arctan model 


0.089 


-0.216 


-0.015 


-0.0298 


-0.0036 



TABLE I: Numerical values of spectral indices, their running 
and the tensor-to-scalar ratio for the three models considered. 
All quantities are evaluated at k, = 0.01/iMpc -1 . 

of those cosmological parameters that do not describe the 
initial perturbations, one would like to know whether 
robust results can be obtained using simple forms for 
the initial power spectra rather than detailed inflationary 
predictions. Therefore, in the context envisaged in this 
section, the observer does not use the assumption that 
inflation is the correct underlying theory, other than to 
motivate the restriction of the scalar perturbations to be 
adiabatic. 

It is common practice to assume a power-law shape 
for the spectrum specified by an amplitude and a spec- 
tral index. The reasoning for this parameterization is its 
simplicity. In the absence of any physical model for the 
generation of fluctuations, one assumes that there is no 
distinguished physical scale in the primordial power spec- 
tra. In order to allow for mildly scale-dependent power 
spectra, a running of the spectral indices can be included. 
This leads to the following shape 

/ l \ «flt+3"flt ln(fe/fe») 
Vnt(k)=V ti t(h)(—J , (10) 

where ngt is either rtg — 1 or n^. 3 The pivot scale fc„ 
is the scale at which all the quantities are evaluated. A 
useful way of viewing Eq. I|1(J|I is that it is the first terms 
of a Taylor expansion of m'P(fc) in lnfc about the pivot 
scale, which draws one's attention to the possibility of 
using other expansions. 

The simplest assumption would be to take both spec- 
tra as constant (scale- invariant) . Models with ns — 1 = 
and tit = do in fact provide acceptable fits to recent 
CMB data (for sufficiently low or zero tensor amplitude) ; 
thus, if we decide to ignore inflation for the moment, 
there is no reason from CMB observations alone to in- 
clude a tilt. Ref. quotes n s - 1 = -0.07t ;^ at 95% 
confidence level, and the addition of large-scale struc- 
ture data greatly tightens the constraint without altering 
the conclusion that ng = 1 is allowed. Current observa- 
tional constraints on the running of the spectral index 
are far weaker than the magnitude predicted in popu- 
lar inflationary models. Thus far, only upper limits on 
the contribution of gravitational waves have been derived 
Il5l ITf^l . Some of these limits suffer from the problems 



no evidence that they are required. 
2 A general analysis would also have to consider vector modes and 3 Note from the definition of the spectral index that ns(k) - 1 ^ 
the various possible isocurvature modes, but at present there is rifl t + ^Qflt ha(k/k») away from the pivot scale. 
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described below Eq. and use strong priors on some 
of the other cosmological parameters. Translating the re- 
sult of Fig. 5 of Ref. 7] (r < 0.5 at 95% confidence level) 
to our notation gives R = 9r/25 < 0.2, while Ref. [lg 
gives a weaker constraint also consistent with R = 0. Let 
us also remark that the majority of recent papers esti- 
mating parameters from the microwave background have 
done so under the assumption that the scalar spectrum 
has a power-law shape and that there is no contribution 
from tensor perturbations (R = 0) 

The question of how far power spectra expansions 
should be taken, and how accurately their coefficients 
need to be computed, obviously depends on the accu- 
racy and dynamic range of observations. For present ob- 
servations an accuracy level of ten percent or better is 
certainly required. Ultimately, Planck will measure mul- 
tipole moments Ce from I of 2 to about 2000, correspond- 
ing to Alnfc/fc* ~ 3.5 on either side of a central pivot 
fe* . It is rather unclear how accurately the multipole mo- 
ments need to be represented at the extremes (cosmic 
variance intervening on large scales and the damping tail 
removing the signal on short scales), but in the centre an 
accuracy of better than one percent is certainly desired 
(see e.g. Ref. ^3)- 4 ^ one then further assumes that 
Planck data will be combined with high-accuracy galaxy 
correlation data, the fc-range might extend to around 
Alnfc/A;* ~ 6 (corresponding to fc max — 30/i Mpc -1 ), 
though the nonlincarly-cvolved galaxy power spectrum 
on short scales is unlikely to be amenable to extremely 
accurate multi-parameter estimation. The choice of pivot 
scale fc* is important as the difference between the fitted 
and the true power spectrum produces an error that runs 
as we move away from the pivot scale. While a careful 
tracking of error covariances should lead to results inde- 
pendent of the choice of pivot, those covariances should 
be minimized to a good approximation by aligning fc* 
with £*, the multipole where we expect the observational 
errors to be least. One can use the approximation 0] 



Ho 



2 1 + 0.084 In ft n 



n A = i, (ii) 



where H = h/3000 Mpc -1 to carry out this alignment. 

Having described what are the typical errors in the 
multipole moments, Error (Cg), we now need to link 
this quantity to the error in the power spectrum itself, 
Error("P), since this is the quantity calculated in prac- 
tice. We assume throughout this paper that an error 
in our determination of the power spectrum propagates 
directly to an error in our determination of the C/s since 



Q = 4tt / dlnkV(k)[A e (k)}' 



(12) 



4 We note that current implementations of CMBFAST ^3 an d CAMB 
[Til have a target accuracy of one percent, so there is presently 
nothing to gain by demanding power spectrum accuracy much 
higher than this. 



where A^(fc) is the £-th momentum of the temperature 
fluctuations. In other words, we assume that Error(C^)~ 
Error 

Another question is how an error in the power spec- 
trum propagates to an error in the estimation of the 
cosmological parameters. In general the Fisher matrix 
formulation is needed to estimate how well a given ex- 
periment can measure the parameters; the error in the 
cosmological parameters is not simply related to the er- 
ror in the power spectrum as there are many parameters 
and lots of degeneracies amongst them. The requirement 
Error(7 : ')~ 1% for Planck is a very stringent condition. 
In particular, it does not imply that parameter estimates 
would go astray if we drifted outside our power spectrum 
accuracy criterion; we would expect parameter estimates 
to stabilize some way before the fitted power spectrum 
was within our 1% accuracy everywhere. Our criterion is 
a sufficient and conservative condition for establishing a 
safe procedure: as long as the power spectrum accuracy is 
below 1% everywhere, we are confident that the system- 
atic errors coming from an inaccurate parametrization 
of the initial conditions will not play a role in the data 
analysis of an experiment like Planck. 



B. Accuracy of the parametrized spectra 

We now investigate the systematic errors which might 
arise from assuming that the spectra have perfect power- 
law shapes or, in a more sophisticated version, a constant 
value for the running of the spectral index. The first step 
is to fix the numerical values of the coefficients Pfit^*), 
rifit and ag t . We have no means to calculate them the- 
oretically in the present context. In practice, observers 
determine these coefficients by carrying out a fit to the 
data. Here, we carry out a least-squares fit of Vat(k) to 
"Pnum(fc) to obtain best-fit scale-invariant, power-law and 
power-law plus running spectra. This means that the co- 
efficients 'Pfit(fc*), "-fit and Qifit are those for which the 
quantity 



E 



(13) 



is minimized. We took the fcj to be equally spaced in 
dlnfc and given equal weight. This idealized fitting ap- 
proach will tend to sacrifice accuracy in the centre of de- 
sired range in favour of accuracy at the extremes. Here 
the idea is to test whether in principle the shape of Pat 
can reproduce the true power spectrum over a reasonable 
range in k. This obviously becomes important if, for ex- 
ample, we try to use a power-law shape to fit to a model 
with significant running of the spectral index. The result 
of the minimization procedure for the three models in- 
troduced above is summarized in Table ITT1 These values 
should be compared with the exact ones of Table HJ 

For the first two examples we conclude that the se- 
quence of fitting a constant amplitude, a power-law, and 
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Fitted values 


-*^fit / num 


R 


ns — 1 


TIT 


as 


QT 


Chaotic 


1.05 


0.279 












1.00 


0.285 


-0.054 


-0.036 








1.00 


0.285 


-0.055 


-0.037 


-0.0010 


-0.0007 


False vacuum 


0.98 


0.053 












1.01 


0.051 


0.054 


-0.006 








1.00 


0.051 


0.054 


-0.006 


0.0017 


0.0004 


Arctan model 


1.23 


0.072 












0.95 


0.092 


-0.178 


-0.016 








0.99 


0.090 


-0.238 


-0.020 


-0.0303 


-0.0044 



TABLE II: The ratio of the fitted amplitude to the numerical 
amplitude of the scalars, -Afit/nuni) and the best-fit values of 
the tensor-to-scalar ratio, the spectral indices and their run- 
ning for the three models considered at fc* = 0.01/iMpc - . 
For each model, we present the results for a scale-invariant, 
power-law and power-law with running spectral shape in three 
rows respectively. 



finally a power-law with running provides best-fit val- 
ues which reproduce the numerical values with sufficient 
accuracy. From the observational point of view this is 
reflected in the fact that the best-fit value of R is the 
same in the second and third row and that the fit values 
of the spectral indices are almost the same as well. Such 
a behavior is the experimental evidence that the input 
does make sense. The situation is different for our third 
example, the arctan model. Although there is slow con- 
vergence in the best- fit values of R, no sign of convergence 
can be detected by inspection of Table ITT1 in the spectral 
indices. This is confirmed by a comparison with the nu- 
merical values of Tabled e.g., the fitted spectral index of 
the tensors is less precise in the third row than in the sec- 
ond, the scalar spectral index is underestimated by 0.036 
by the power-law fit and overestimated by 0.022 including 
running. From the point of view of inflationary param- 
eters, see below, these are large fluctuations. Thus an 
observer in possession of sufficiently-accurate data some- 
time in the future should conclude for our third example 
that more parameters have to be introduced in the fit, 
before any physical meaning can be extracted from the 
best- fit values of the spectral indices. 

A large error in the fitted values of the amplitude and 
of the spectral index, even if we are not interested in the 
physics of inflation for the moment, is undesirable for two 
reasons. Firstly the overall amplitude of scalar perturba- 
tions is a quantity that we hope to measure from MAP 
and Planck at the percent level. Thus we would prefer 
to relate it to a physical quantity, namely the amplitude 
of super horizon density fluctuations. The second reason 
comes from considerations of large-scale structure data, 
where it is customary to include a linear bias parameter, 
6, to account for the overall normalization of the matter 
power spectrum. If we simultaneously fit to the CMB, 
then we can only assign any physical meaning to b if we 
are certain that the amplitude of scalar perturbations is 
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FIG. 2: Error curves for various fits to the scalar power 
spectrum for the false vacuum model. While the power-law 
fit is acceptable for fitting to present data, neglecting running 
affects the estimate of the power spectrum amplitude at the 
pivot point at the percent level. 

correct. In addition, an inaccurate estimate of the am- 
plitude and the tilt could spoil a consistency check of 
structure formation based on measurements of ag. 

Having determined the coefficients, the second step is 
now to compute the error. We define this by 

Error(7>) = ( - 1 ) X 100% . (14) 

\ / num / 

In the following, we give three examples. 

In Fig. |21 we plot the error in the scalar power spectrum 
in the case of the false vacuum model. The best-fit scale- 
invariant spectrum is a poor fit for this particular model, 
while the best-fit power-law spectrum improves things 
greatly, keeping the errors below 2% which is more than 
adequate for present CMB data and marginally adequate 
for Planck. The large effect of the tilt is due to a long 
lever arm in wavenumbers p(il | ; the error being of the or- 
der (ng — 1) In fc /fc* , even a small tilt can have a significant 
effect if the data span several decades in wavenumbers. 
We can further see that with the inclusion of running, 
Pat now reproduces the power spectrum in great detail. 
This is because the correction to the spectrum is of or- 
der as In 2 k/k* which, for the running of this example of 
as ~ 0.002, gives a significant effect though the correc- 
tion is much smaller than that from the tilt. 

Given a set of observations, the importance of running 
is tested by including it in the fit and examining whether 
the fit improves significantly. In the absence of any the- 
oretical prejudice, one might well hope to detect signifi- 
cant running at high significance. However, some of the 
simplest inflation models predict running of at least an 
order of magnitude below what even Planck can achieve 
[2l| . In that case there will be no significant detection of 
running, and marginalizing over the running permitted 
by the observations may lead to a significant inflating of 
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FIG. 3: As Fig. but for the arctan model. For fitting to 
data of the present quality, the inclusion of the running is 
required. 



errors on other parameters. While combining short-scale 
observations with the microwave background may give a 
stronger lever-arm in constraining running, this may well 
turn out to be a parameter for which it is desirable to in- 
vestigate imposing a strict theoretically-motivated prior 
to compare with a free fit. Further, even if running is de- 
tected at high significance this problem then resurfaces 
concerning the running-of-running. 

As we have already concluded from the discussion of 
the best-fit values in Table [HJ there exist models for 
which a power-law fit to the spectrum does not provide 
a good description in contrast to the above example. In 
Fig. |3| the error in the scalar power spectrum for the arc- 
tan model is displayed. Including the running is neces- 
sary for the present accuracy of CMB experiments. Now 
the effect of running is comparable to that of the tilt. 
We see that more parameters (e.g. running of the run- 
ning) would be necessary to reproduce the power spec- 
trum with 1% accuracy. We actually have to add more 
and more parameters until the spectrum starts to con- 
verge (see also the discussion of Table |Hj| , or consider 
using a different spectral shape. 

For the tensors in the case of the chaotic model, we see 
in Fig.^that the spectrum is poorly fitted by the scale- 
invariant spectrum. However, the accuracy requirements 
on the tensor spectrum are less stringent — the tensor 
amplitude is generally less that the scalar amplitude and 
so the required absolute error in Th is also less [~ R x 
Error^^)]. Thus, a typical inflationary tensor spectrum 
is well described even by a scale-invariant spectrum for 
present-day experiments, though there is no reason not 
to describe it with the same sophistication as the scalar 
spectrum. For future CMB measurements, the inclusion 
of a tilt is sufficient in this example. 

To answer the main question of this section, we can 
expect to obtain robust estimates for the cosmological 
parameters for a restricted class of inflationary models 
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FIG. 4: As Fig. |3but for the chaotic inflation model tensor 
spectrum. Although the percentage error is large for the scale- 
invariant fit, the absolute error is small compared to the scalar 
spectrum, and so the scale-invariant fit is still acceptable. 



using the fitting procedure described above. However, 
there exist models where this is no longer true. In the 
following sections, we specify the criteria which define 
this class of models. 



IV. PREDICTIONS OF SLOW-ROLL 
INFLATION 

In this section, we restrict our considerations to the 
class of slow-roll models of inflation. The advantage is 
that we can now predict the shape of the power spec- 
tra and link the parameters characterizing these spectra 
to the physics of inflation. There has recently been re- 
newed progress in the accurate calculation of inflationary 
perturbations by analytical techniques, including a com- 
putation of the power spectra to arbitrary order in the 
slow-roll expansion for single-field inflation by Stewart 
and Gong |22|. and a computation at higher order for 
models that may violate one of the slow-roll conditions 
p3l |24| . We utilize the Stewart-Gong results here as 
they have the most general applicability, extending them 
with an explicit evaluation of higher-order terms for the 
tensor spectrum. 

The background evolution can be described in terms 
of the horizon-flow parameters {e n } |3- Starting from 
e = H(Ni)/H(N), where l/H is the Hubble distance 
and N = ln(a/ai) the number of e-folds since some initial 
time t[, the set {e„} is defined by 



These parameters can be easily related to various def- 
initions of the slow-roll parameters. Setting n = 1 we 
find t\ = — dln-ff/dlna, which is nothing but the slow- 
roll parameter e of Refs. [25L |2(j. The parameter r\ of 
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Refs. |25|, |26(, which is usually defined to measure the 
deceleration of the inflaton field, enters as e 2 — 2e — 2rj. 
The third slow-roll parameter, £, is contained in 6263 = 
4e 2 — 6er] + 2£ 2 . In this notation, all the e n are typically 
of the same order of magnitude. Inflation takes place 
provided ei < 1. Slow- roll inflation is defined by the 
condition \e n \ <C 1, for all n > 0. 

A measurement of the horizon-flow parameters, at a 
specific moment during inflation, would immediately pro- 
vide us with a value for the inflaton potential V and its 
derivatives with respect to the inflaton field <fi (denoted 
by a prime in what follows) for any single-field inflation 
model. For example, from H and £i,£ 2 , and £3 we can 
calculate the potential and its first two derivatives ex- 
actly, 



V 

V 

V" 
3H 2 



8tt 
3m Pl H 2 
' (4tt) 1 /2 



2ei 



£2 
2 



1/2 



2e 2 



_ £1 
3 

5ei£ 2 
6 



£2 
6 

i 

12 



£2f3 

6 



(16) 
(17) 

(18) 



If £3 cannot be determined and the horizon-flow param- 
eters are small compared to unity, we can still estimate 
V" by keeping the leading terms only. 

For slow-roll models we can invert this procedure and 
estimate the horizon-flow parameters. At leading order 
in these parameters we find: 



H 



£1 



(-2 



£2 £3 



8tt 



" 6 P1 
16tt 

47T 

4 

32tt 2 
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v) ' 
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~V"'V' 
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V 



(19) 
(20) 
(21) 

.(22) 



To give an example, for chaotic inflation with the poten- 
tial V oc 7 we find ei ~ j/AAN and e 2 ~ e 3 ~ 1/AiV, 
where AN denotes the number of e-folds before infla- 
tion ends. Chaotic inflation is a simple model where 
the higher horizon-flow parameters are of the same or- 
der of magnitude as lower ones. In the case of power- 
law inflation (a oc t p ) where the potential is given by 
V oc exp[— / p) 1 / 2 4> / mp\\, we recover the exact result 
€\ = 1/p and £2 = £3 = 0. 

The power spectra of scalar and tensor perturbations 
can be obtained approximately using analytic techniques. 
One expands the power spectra about some particular 
wavenumber fc* , and then computes the coefficients using 
the slow-roll expansion or some other scheme of approxi- 
mation. This amounts to a double approximation. Given 
that we need to cover several orders of magnitude in fc, 



the most appropriate expansion variable is In k, giving 

P(k) , / k\ 02, a / k 

-^o + axlnf- +-ln 2 - 



(23) 



The next step is to establish an expression for the coeffi- 
cients a n , which can be obtained either with help of the 
slow-roll expansion |E El IH IH IH Hlor the meth- 
ods of approximation developed in Refs. |23l 0] • Since 
the former covers a more general class of inflation mod- 
els then the latter, we focus on slow-roll inflation in the 
following. We will use the term first-order to refer to re- 
sults including all terms up to order e m and second-order 
if one goes to terms including e 2 m . 

The normalization of the power spectra is set by the 
expansion rate during inflation, H, and the parameter 
ei, namely 

Vko(K) - — (24) 



7reimpj 
16H 2 



(25) 



'pi 



where H and ei are evaluated when aH = fc, during 
inflation. The scalar amplitude has been calculated up 
to first-order in the slow-roll parameters by Stewart and 
Lyth [23, and recently up to second-order by Stewart 
and Gong 22]. These calculations are sufficient to al- 
low calculation of an infinite, though incomplete, set of 
expansion coefficients of which the first few are given by 



aso 



1 



-2(C 

'c 1 - c 



l)ei-Cc 2 



7tt 7 

12 ' 



2C 2 + 2C 



f -5 



£l£2 



£ 



-iC 2 



asi 



as2 



-2ei 
-CS, 



4e 2 



£ 2 £ 3 , (26) 
e 2 + 2(2C+l)e 2 + (2C-l)eie 2 
2 - C£2£3, (27) 
2£i£ 2 + £ 2 -£ 2 £3, (28) 
—0.7296. For the tensors, the 



where C = + In 2 — 2 
corresponding set is as follows: 

= l-2(C + l)ei + (2C 2 
+ (-C 2 - 2C 



OTO 



2C + \ - 5 



12 



-2 )e 1 e 2 , 



OTl 

ax2 



4, 2 



(29) 

2e 1 + 2(2C*+l)e 2 -2(C + l)e 1 e 2 , (30) 

2£i£ 2 . (31) 

We have presented for the first time the C(e 2 ) terms in 
the tensor amplitude which we obtained along the lines 
ofRef. [12. 

The coefficients a n for n > can also be obtained by 
successive differentiation of the first term of the expan- 
sion 

d n [V(k)/V (K)] 



dln n fc 



(32) 



k=k, 



1 



1 - ei dN 



Vo(h)ao(K), (33) 
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Slow-roll values R ns — 1 tit 



as 



CKT 



10 



Chaotic 


0.285 


-0.055 


-0.037 


-0.0010 


-0.0007 


False vacuum 


0.051 


0.054 


-0.006 


0.0017 


0.0004 


Arctan model 


0.089 


-0.221 


-0.014 


-0.0291 


-0.0041 



TABLE III: Slow-roll values of spectral indices, their running 
and the tensor-to-scalar ratio for the three models considered. 
All quantities are evaluated at k — 0.01/iMpc _1 . 



where we used the "horizon crossing" condition fc* = k = 
aH to obtain the second line. From Eqs. (|15f) and i|33|) 
we see that the leading contribution to a n is of order 
ej^ (where means any terms containing n of the e, 
not necessarily all the same). If ao has been written to 
first order, differentiation yields a\ to second order, a 2 to 
third order and so on. Note that the coefficients of the 
Taylor series, Eq. (|23[1 . always feature an increasing num- 
ber of powers of the slow-roll parameters, so in practice 
convergence of the Taylor series is governed by the size 
of e m lnfc/fc*, which in principle needs to be small for all 
values of m > 1. Thus the series can still be strongly con- 
vergent even if ln/c/fc* exceeds one, as it will for typical 
upcoming experiments. 

Let us now calculate the spectral indices and their run- 
ning in the slow-roll approximation up to second order. 
For this purpose, it useful to calculate the logarithm of 
the power spectrum 



hi 



_ = 6o + 6l In^-j +T In 



(34) 



Exponentiation of Eq. I|34|) automatically enforces the 
positive definiteness of V(k) and allows us to directly 
link the first coefficients b n to the spectral indices and 
the runnings, because 

&Si = ns ~ 1, b T i = n T , b S 2 = a s , b T 2 = a T - (35) 
The equivalent expressions to Eqs. I|2tj|) - (|31|l are 

b so = -2(C+l) ei -Ce 2 + (-2C+^-7)e? 

-C 2 - 3C + 3f - 7) e l£2 + - l) 4 



(-IC 2 + £) e 2 e 3 , 



(36) 



-2ei - e 2 - 2et - (2C + 3)eie 2 - Ce 2 e 3 , (37) 



bs2 = — 2eie 2 — e 2 £3, 
for the scalars, and 

b TQ = -2(C + l)e 1 + (_2C+£-7)<7 

+ (-C 2 -2C+f^-2) £l e 2 , 

frn = -2ci - 2e? - 2(C + l)eie 2 , 
&T2 = — 2eie 2 , 

for the tensors. 



(38) 



(39) 

(40) 
(41) 
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FIG. 5: Fitting the slow-roll shape to the arctan model. The 
errors should be compared with the errors in Fig. |3 For this 
model, the second-order slow-roll shape provides a better fit 
that the power-law plus running shape. 



Finally, the ratio of amplitudes of scalars and tensors 
at the pivot point is 



R 



16ei 



1 + Ce 2 + (C 



5 eie 2 



+ 



1 



94" £ 2£3 



.(42) 



This becomes the well-known "consistency condition of 
inflation" R = — 8nx at leading order, which holds for 
single-inflaton- field slow-roll models. The values of the 
ratio R, the spectral indices and their running, computed 
in the slow-roll approximation for the three models en- 
visaged in this article, are summarized in Table ITTT1 The 
values of the horizon-flow parameters were obtained nu- 
merically, though an actual reconstruction may also fea- 
ture a slow-roll approximation in relating those to the 
inflationary potential. 



V. DOES THE SHAPE OF THE FITTED 
SPECTRA MATTER? 

In the preceding section, we have shown that the shape 
of the slow-roll spectra does not coincide with the shape 
of the fit of Section III. From a theoretical point of view, 
it is clear that the former should be used not only to 
predict the spectra but also to fit real data. For many 
choices of parameters the difference between the shapes 
is not significant, but there are also models where this 
difference can be important. 

An example is given in Fig. [3 where we plot 



Error(7 ? , 



if it 



1 x 100%, 



(43) 



for the arctan model of Section II. In this equation, P^r 
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Slow-roll fit 


•Afit/nun 


R 


ns — 1 


TIT 


as 




Chaotic 


1.05 


0.279 












1.01 


0.283 


-0.056 


-0.037 








1.00 


0.285 


-0.055 


-0.037 


-0.0010 


-0.007 


False vacuum 


0.98 


0.053 












1.01 


0.051 


0.051 


-0.006 








1.00 


0.051 


0.055 


-0.006 


0.0017 


0.0004 


Arctan model 


1.23 


0.072 












1.07 


0.082 


-0.210 


-0.016 








1.00 


0.089 


-0.213 


-0.019 


-0.0289 


-0.0044 



TABLE IV: As in Table E] but for the spectral shape that 
is predicted by slow-roll inflation. 



is found by considering 



<< nlu( Pl+f ln 2 (A), ,44, 



and calculating the three coefficients Co, c\ and ci by 
minimizing the quantity 



E 



(45) 



Comparing the slow-roll fit of Fig. [5] with the power-law 
fit of Fig. |21 we can see that the slow-roll shape does 
indeed provide a better fit in this case, keeping the er- 
ror below 1% for most of the range. Thus the power 
spectrum shape can make a difference, and there exist 
models where fitting with the power-law instead of the 
slow-roll shape can lead to significant errors ( defined by 
the criterion of Sec. Ill- A ). 

However, one cannot conclude that the slow-roll shape 
necessarily gives a better fit in general. An example 
where the slow-roll fit converges slower than the power- 
law fit is the chaotic model, although the difference is 
not significant in that case. For power-law inflation the 
slow-roll shape will actually fare less well. 

A second step is to go from the coefficients Co, c\ and ci 
to the characteristic parameters of the primordial spec- 
tra. This can be done by means of the relations 



(ns-l)Sf = - + 0(e 3 n ), 
co 

(«s)2? = --4 + 0(4), 

co cf) 



(46) 
(47) 



and analogous equations for the tensors. The coefficient 
R can be obtained as 



ofit = cot 
cos 



(48) 



The results are summarized in Table IIVI This table 
should be compared with Tables III and HT1 Fitting a dif- 
ferent shape has now the effect that the parameters of 



the arctan model converge, in contrast to the power-law 
fit. 

Fitting the coefficients c n allows us to test the con- 
sistency relation of inflation, and thereafter constraining 
cit and c 2 t according to Eqs. lj50"ll and allows us to 
measure the inflationary parameters. 

Having shown that there exist situations where the 
shape matters, we wish to find the region of the param- 
eter space in which the difference between a power-law 
shape with running and the shape predicted by slow-roll 
inflation is significant. For this purpose, we define the 
estimator 



Vsr - Vat 

(V S r + Pfit) /2 
2 



x 100% : 



n 

~2" ! 



(49) 



y jln^-) x!00%, (50) 



"~ x i ■ 3 I k 



where n stands in for ns — 1 or nx- Note that this estima- 
tor presumes that the two fits generate the same values 
for the amplitude, spectral index and running, whereas 
in practice a different choice of shape will lead to differ- 
ent values. This estimator therefore underestimates the 
differences between the two fits close to the pivot point 
and overestimates them far away from the pivot. 

In Fig. |f)| we plot the contours of the maximum of 
|cr(fc)| in the interval —1.5 < log 10 (fc/fc*) < 1.5 in 
the (ns — l,as) plane. Its shape can be understood 
most easily from the approximation Eq. (|50|l . We con- 
clude that within the ranges ns — 1 G [—0.05, 0.05] and 
as G [—0.015,0.015], shape should not matter even at 
the accuracy level of Planck. For present CMB experi- 
ments this plot suggests that as long as |ng — 1| is within 
the range shown in Fig. H3 shape is not an issue if the 
running is at most of order 0.01, which is the case for 
a wide class of inflationary models (similar constraints 
should be assumed to hold true for higher corrections as 
well). 

A significant difference between the two fits at a given 
observational accuracy is a clear indicator that higher- 
order terms may be important, as it is those which give 
the difference between the two expansions. To be certain 
of robust results, an attempt should be made to estimate 
these higher-order terms, either by extending one or both 
expansions to seek convergence between them or by re- 
sorting to fully numerical analysis techniques. 



VI. ACCURACY OF SLOW-ROLL ANALYTIC 
SPECTRA 

In the previous section we showed that spectral shape 
can matter and therefore that it is important to take the 
predictions of slow-roll inflation into account if we are 
interested in the physics of inflation itself. Before dis- 
cussing how to extract the inflationary parameters we 
study the accuracy of the slow-roll approximation at sec- 
ond order. First studies of the accuracy of the slow-roll 
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n 



FIG. 6: The region of fitted spectral indices and runnings 
in which the difference between the power-law shape and the 
slow-roll shape, estimated by \cr\, is within 1% and within 
10%. 



expansion can be found for the amplitudes in Ref. |j| 
by comparing to numerical results, while in Ref. |2iJ the 
first-order expressions for the amplitudes and the spec- 
tral indices has been tested by comparison to analyti- 
cal results for power-law inflation. Here we extend these 
studies to the full power spectrum at second order. We 
define the error of the slow-roll power spectrum as 

Error(^) = ( — ^ - 1 ) x 100%, (51) 

\ f num / 

where V sr is given by Eqs. J2U) and P5 Jh J5T |) . In these 
expressions the values of H, e\, €2 and e3 are computed 
numerically for the three models of Sec. II. 

Looking at the chaotic inflation model first, we can 
see from Fig. that the error curves resulting from slow- 
roll predictions generally have the property that they are 
most accurate close to the pivot point (in terms of ampli- 
tude and spectral index) and that the error increases as 
we move away from the pivot point. We can also see that 
the second-order expressions can improve the accuracy of 
both the scalar and tensor power spectra to within Planck 
requirements, whereas the accuracy of the corresponding 
first-order expression would be at best marginal. This 
improvement is mostly brought about by the inclusion of 
the running. 

The tensor spectrum of Fig. is determined more ac- 
curately than the scalars. We have observed that this is 
typically the case. Since the accuracy requirement upon 
the tensors is less than on the scalars, it is the scalars 
upon which attention should be focused. 

Next we turn to the false vacuum inflation model. Note 
immediately from Fig.[8]that the second-order expression 
improves both the shape of the power spectrum and the 
accuracy of the amplitude at the pivot point itself. The 
first-order expression is good enough for present experi- 
ments in this example, but not for MAP and Planck. 
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FIG. 7: Scalar and tensor error curves for the chaotic inflation 
potential. The pivot scale crosses the Hubble horizon 55 e- 
folds before the end of inflation. We see an improvement in 
accuracy from the first to the second-order expressions. The 
tensors have better overall accuracy than the scalars. 



Finally, for the arctan model we see in Fig. [U] that al- 
though ei is small and 62 and 63 are still in agreement 
with the slow-roll conditions, the effect of the second- 
order correction is very important. The first-order ex- 
pression is not sufficient for MAP. In this example, the 
first-order expression also produces a significant error in 
the amplitude at the pivot point. For Planck the plot 
suggests that the third order is necessary. 

It is of course impossible to study the accuracy of all 
possible models of inflation in this way. We therefore 
need a more general estimator for the accuracy of the 
slow- roll expansion in the parameter space e„. The differ- 
ence between the slow-roll expansions of V(k) and lnP(fc) 
is such an estimator. We define the error at a given order 
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FIG. 8: Scalar error curve for the false vacuum inflation 
model. Again, we see an improvement in accuracy from the 
first to second-order expressions which helps to correct the 
amplitude at the pivot point. 



FIG. 9: Scalar error curve for the arctan potential. We see 
an improvement from the first to second order as well as a 
significant correction to the overall amplitude at the pivot 
scale. 



n to be 





- exp 










+ exp 









xl00%, 



where the coefficients and hi are taken at order e^. 
The interpretation of this expression is that it gives the 
smallest fractional amount by which the worse of the 
two expansions departs from the true power spectrum, 
namely half the distance between the two estimates. This 
interpretation justifies the absence of a factor 1/2 at the 
denominator in Eq. (|52|) . 

This expression is of order e^ 1 and therefore is an in- 
dicator of the importance of orders that have not been 
included. Moreover it has the same typical behaviour of 
the errors as one goes away from the pivot point, and 
we also find that it estimates the orders of the errors 
for the examples of Sec. II correctly. We expect that 
this estimate typically works well although there exists 
the possibility of fine-tuning models such that a n is not 
a good estimator. In the following we study the maxi- 
mum of the error in a suitable interval of wavenumbers, 
because a large error in a small range may spoil an oth- 
erwise accurate fit. We therefore maximize crx(k, £1,62) 
and <J2(k, £1, e 2 , e 3 ) over —1.5 < log 10 (fc/fc») < 1.5. This 
is certainly conservative but is a good indicator of when 
robust results are expected. 

The upper panel of Fig. I1UI shows the error in the t\- 
£2 plane, maximizing over —0.1 < £3 < 0.1 (the arctan 
model actually lies outside this range). The scalar error 
contours are elongated along the direction £1 = —£2/2, 
which corresponds to ns = 1 at first order. In the top 
left corner o\ becomes independent of the dominant con- 
tribution proportional to ln/c for ns = 1. For 172 there is 
a similar cancellation of the In 2 k contribution for models 
close to ns = 1, which explains the shape of the contours. 



These elongated shapes are therefore a feature of our es- 
timator <r n ; they do not reflect a proper estimate of the 
error in the top left corner as other higher-order terms 
not considered would spoil these cancellations. With the 
exception of that top region, we see that, as expected, 
the second-order expressions extend the area of parame- 
ter space meeting a specified accuracy requirement. 

It is useful to examine these results in the (ns — 1)—R 
plane via the transformation 

n s -l = -2ei-e 2 , (53) 
R = 16£i (54) 

shown in the lower panel of Fig.^| We use the first-order 
relations also for the second-order error contours here; 
the error made by this can be neglected for the present 
purpose. The restriction that we put on £3 gives rise to 
values for running in the range as € [—0.023,0.14] for 
the displayed region of parameter space. The first-order 
expression gives errors within 10% in the region given 
approximately by —0.15 < ns — 1 < 0.1 and R < 1.5. 
The second-order slow-roll expression gives an accuracy 
better than 1% in a somewhat smaller range of parameter 
space (-0.1 < n s - 1 < 0.05, R < 1.0). 

It is important to stress that these regions are very 
conservative as we maximize the error over both £3 and 
wavenumber. The conclusions of small errors in parame- 
ter space regions is therefore very robust, and indeed the 
errors are likely to be within acceptable levels even for 
many models lying outside our contours. 

An important limit is when £1 is very small, since a 
broad class of inflation models belong to this category, 
e.g. false vacuum dominated inflation gives rise to tiny 
e\. When £1 < 0.001, then the tensor spectrum will have 
no effect on the low-£ portion of the Ct curves at the 1% 
level, see Eq. Q42[l. At this point the tensor C/s drop 
out of reach and we can no longer measure H during in- 
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FIG. 10: These panels show the error estimate <ti for the 
slow-roll expansions at first order (thin lines) and 02 at sec- 
ond order (thick lines). The upper panel is as a function 
of horizon-flow parameters, while the lower panel transforms 
this into the (ns — 1)—R plane. 



flation and t\ separately, see Eq. (|24fl . The scalar power 
spectrum, Eqs. l(2l)|> - J2HJ), now reduces to a function of 
'Piza(k*), £2 and £263, where the last two parameters de- 
termine ng — 1 = —62 — Ce2C3 and as = —£2^3. In Fig. 1111 
we plot the error of the second-order power spectrum, 172, 
in the (ns — 1) " <^s plane. The transformation between 
the (ns — 1) - as plane and the £2 - £3 plane is nonlinear 
and singular at £2 = for any £3. All corresponding mod- 
els have ns — 1 — a s — 0. Moreover, in the vicinity of the 
line as = (ns — 1)/C the value of £2 becomes arbitrarily 
small, and thus £3 can be huge. Therefore, in the vicinity 
of the dashed line the estimator cr 2 is misleading, because 
it gives a small error even for models which violate the 
slow-roll condition £3 <§; 1. Nevertheless, the conclusion 
is that fairly weak running < 0.02 can be accurately (1%) 
described by a slow-roll expansion with tiny ei. 



FIG. 11: The error estimate 02 in the (ns — 1) - as plane, 
with ei <C 0.001. The dashed line is as = (ns — 1)/C, in the 
vicinity of which the error estimate can be misleading. 



VII. TESTING SLOW-ROLL INFLATION 

We end with a proposal of how to proceed with 
testing slow-roll single-field inflation using future high- 
accuracy data. The corresponding data analysis pipeline 
is sketched in Fig. ^| The inputs are the CMB data 
and a cosmological model (e.g. ACDM). The first step 
should be to determine the cosmological parameters un- 
der the assumption that the power spectra of scalar and 
tensor perturbations are given by a power-law with run- 
ning of the spectral index, see Eq. I|10[) . One should check 
the convergence of the values of all cosmological parame- 
ters as one fits scale-invariant, power-law, and power-law 
with running spectra, as discussed in Sec. lIIII One should 
continue to refine the power spectrum shape (adding in 
running of running etc) until the new power spectrum 
parameter is found to be consistent with zero. At this 
point one has the choice to neglect this final parameter, 
and this seems a sensible option. We call the order of 
this truncated power spectrum O p \. In a similar manner 
one should also check the convergence of the cosmolog- 
ical parameter estimates while fitting to the data using 
scale-invariant, first-order and then second-order slow- 
roll shapes, up to order O sy . 

One should find \O p \ — O sr | < 1, with O v \ — O SI being 
the most likely case. If we also find consistent estimates 
of the cosmological parameters then clearly the choice of 
power spectrum shape doesn't matter. If O p \ ^ O SI but 
the cosmological parameter estimates are convergent and 
consistent with each other, then we have some evidence 
that a particular power spectrum shape may be preferred. 
Figure might be used to check whether the extracted 
spectral indices and runnings are expected to give rise to 
a significant difference between the two fits. 

If there is no convergence using one or both of the 
power spectrum shapes, or if the different power spec- 
trum shapes lead to significantly different estimates of 
the cosmological parameters, then there is either a sig- 
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FIG. 12: Suggested pipeline to test slow-roll inflation and 
estimate its parameters. 



nificant problem in the assumed cosmological model or 
the shape of the spectrum is completely different from a 
power-law, e.g . a pronounced bump or a step at a privi- 
leged scale 29]. Presuming the latter, within the context 
of single-field inflation, the optimal strategy is a direct 
estimation of the inflationary potential from the data it- 
self, without using intermediate approximations such as 
the slow-roll expansion, as described by Grivell and Lid- 
die [3(]]. 5 Such a calculation must simultaneously fit all 
parameters, and so will also test whether the results are 



5 The inflationary potential is parametrized, for example by a Tay- 
lor series, and the scalar and tensor power spectra are obtained 
by solving the m ode equations and fed into a Boltzmann code 
such as CMBFAST or CAMB ^|. The only approximation is 
the validity of linear perturbation theory. The result is an un- 
biased estimation of the inflationary potential with automatic 
generation of the error covariances of the potential param eter s 
amongst themselves and with the cosmological parameters 
Other considerations of single-field inflation beyond slow-roll are 



consistent with a flat universe; the simplest models of 
inflation predict fitotal = 1 ± 10~ 5 , though realistic ex- 
periments will be orders of magnitude larger in uncer- 
tainty. If so the data are consistent with inflation, but 
single- field slow-roll inflation would be ruled out. 

If satisfactory convergence of the cosmological param- 
eters is achieved then the next step is to check whether 
footed is consistent with one. If this test is failed then 
slow-roll inflation is excluded and we need alternative 
physics. If the Universe is consistent with flatness, slow- 
roll inflation can now be taken very seriously. In the 
previous section we have shown that the power of fluc- 
tuations can be predicted at the required level of accu- 
racy in a large region of parameter space favoured by 
present CMB observations. Once slow-roll inflation has 
been adopted as a working hypothesis, fitotal snou ld be 
fixed at unity and not varied in any parameter fits. 

We can now test the consistency relation and then esti- 
mate the inflationary parameters. In principle one could 
use either expansion [Eq. I|1(J|) or Eq. Q44[l] if it has been 
successful, and even if O p \ ^ O SI the inflationary in- 
formation contained within them should be equivalent. 
However presuming it is available it makes best theoreti- 
cal sense to use the slow-roll fit. The approach is sketched 
in the lower tree of the pipeline in Fig. IT51 

The first step is to check whether a tensor contribu- 
tion can be detected at a significant level. If not, then 
there are no means to fully check the specific predictions 
of slow-roll inflation. However, this means that an up- 
per bound on the tensor-to-scalar ratio R is provided by 
the CMB data. Assuming slow-roll inflation we can use 
the consistency relation Eq. I|42[l [or its first-order version 
Eq. I|54[l] to obtain an upper bound on e\. Then we ne- 
glect all ei terms in Eqs. Ij26(l - H31(l . allowing an estimate 
of £2, £3 and the normalization of the scalar power spec- 
trum iJ 2 /7reim P1 . Together with the upper bound on £1 
this gives an upper bound on the scale of inflation H. 
Figure HT1 might be used to estimate the theoretical error 
in the measurement of £2 and £3. If the estimates for | e2 1 
and I £3 1 turn out to be larger than the upper bound on 
£1 we can take these estimates seriously. However, if it 
turns out that one of the higher-order parameters is of 
the same order as the upper bound for £1 we cannot con- 
sistently neglect e±. In this case only a banana-shaped 
region in parameter space of the second-order slow-roll 
expansion can be identified. But a warning is required 
at that point; without a detection of tensors it might be 
impossible to distinguish between single-field slow-roll in- 
flation and other models. 

If there is a significant detection of tensors, the next 
step is to test the consistency equation of slow-roll infla- 
tion Eq. H42|) . If this test is not passed, we have ruled 
out single-field slow-roll inflation. If we find consistency, 
the final step is to measure the scale of inflation H and 



given in Refs. I2II3I . 
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the inflationary parameters ei, €2 and £3. By fitting di- 
rectly for these parameters, rather than the coefficients of 
expansion as above, we are now automatically imposing 
the consistency relations between the scalar and tensor 
spectra. This is also important for measurement of the 
cosmological parameters, as it ensures that the uncer- 
tainties are not overestimated (under the presumption 
that slow-roll inflation is correct) . The slow-roll shape is 
the preferred option for carrying out this final parameter 
determination, and this is also the determination which 
yields the definitive measures of the various cosmologi- 
cal parameters. These might differ from the parameters 
estimated from the power-law plus running fit once the 
consistency conditions are imposed. In particular the 
uncertainties should tighten as the inflationary predic- 
tions are more specific than fitting free power-laws plus 
running. The systematic uncertainty from theory in the 
measurement of inflationary parameters can now be es- 
timated with the help of Fig. I1UI 

Having analytically reconstructed an inflationary po- 
tential, its validity can be checked by evaluating the 
perturbations generated by the potential numerically, 
which will provide a further estimate of the magnitude 
of higher-order corrections. If these prove significant, 
the numerical results could be used to 'tune' the recon- 



structed potential with the aim of removing any biases 
in estimation of other parameters. Ultimately, analytic 
results obtained the way we describe can be compared 
with a direct numerical reconstruction as described in 
Ref. |30|. with the two methods providing invaluable 
cross-checks on each other. 

We have presented a strategy to measure the most im- 
portant quantity in the context of inflationary models, 
the scale of inflation H. It probes the time scale and 
thus the energy scale of new physics, which requires the 
detection of tensor contributions. Sensitivity to gravita- 
tional waves is mainly provided via high-sensitivity po- 
larization measurements, and it is these which may allow 
us to probe the highest energy scales for the first time. 
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